Two-Channel Kondo Lattice: An Incoherent Metal 
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The two-channel Kondo lattice model is examined with a Quantum Monte Carlo simulation in 
the limit of infinite dimensions. We find non-fermi-liquid behavior at low temperatures including 
a finite low-temperature single-particle scattering rate, the lack of a fermi edge and Drude weight. 
However, the low-energy density of electronic states is finite. Thus, we identify this system as an 
incoherent metal. We discuss the relevance of our results for concentrated heavy fermion metals 
with non-Fermi-Liquid behavior. 
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Introduction. The Fermi Liquid theory of Landau has 
provided a remarkably robust paradigm for describing 
the properties of interacting fermion systems such as liq- 
uid 3 He and numerous metals like aluminum and cop- 
per. The key notion of this theory is that the low ly- 
ing excitations of the interacting system possess a 1:1 
map to those of the noninteracting system and hence 
are called "quasiparticles." In the metallic context, one 
may think of the quasiparticles as propagating electron- 
like wave packets with renormalized magnetic moment 
and effective mass reflecting the "molecular field" of the 
surrounding medium. A sharp Fermi surface remains in 
the electron occupancy function which measures the 
number of electrons with a given momentum, and for en- 
ergies lo and temperatures T asymptotically close to the 
Fermi surface the excitations have a decay rate going as 
lu 2 + 7t 2 (/cbT) 2 , which is much smaller than the quasipar- 
ticle energy, which generally translates into a T 2 contri- 
bution to the electrical resistivity p(T). This theory has 
proven useful in describing phase transitions within the 
Fermi liquid, such as superconductivity which is viewed 
as a pairing of Landau quasiparticles in conventional met- 
als such as Al. 

The Fermi liquid paradigm appears now to be break- 
ing down empirically in numerous materials, notably the 
quasi-two-dimensional cuprate superconductors [jj] and a 
number of fully three dimensional heavy fermion alloys 
and compounds In these systems such anomalies as 
a conductivity with linear dependence on u>,T are often 
observed, as are logarithmically divergent linear specific 
heat coefficients. If the quasiparticle paradigm indeed 
breaks down, this may require completely new concepts 
to explain the superconducting phases of these materials. 
While the Luttinger liquid theory provides an elegant 
way to achieve non-Fermi liquid theory in one-dimension 
(with, e.g., no jump discontinuity in n?, and separation 
or unbinding of spin and charge quantum numbers), this 
results from the special point character of the Fermi sur- 
face. Whether the essential spin-charge separation may 
"bootstrap" into higher dimensions remains unclear 



Among the remaining theories to explain experiment are 
those based upon proximity to a zero temperature quan- 
tum critical point 0] , those based upon disorder induced 
distributions of Kondo scales in local moment systems 
B, and those which hope to explain the physics from 
impurity to lattice crossover effects in the multi-channel 
Kondo model [||. Notably lacking for dimensions higher 
than one are rigorous solutions to microscopic models 
which display non- Fermi liquid behavior (7). 

In this Letter, we present the first rigorous solution 
of the two-channel Kondo lattice model in infinite spa- 
tial dimensions. We find that the paramagnetic phase of 
this model is an "incoherent metal" with finite density of 
states at the Fermi energy and finite residual resistivity. 
The excitation spectrum is non-fcrmi liquid like; in par- 
ticular, there is no discontinuity in n^, a finite lifetime for 
electrons at the Fermi energy, an ill defined quasiparticle 
mass, and a linear in temperature saturation of the elec- 
trical resistivity. We find that physical quantities may 
be suitably scaled with a lattice Kondo scale To that is 
significantly enhanced over the impurity limit. We dis- 
cuss the possible relevance of these results to understand 
transport properties of concentrated heavy electron ma- 
terials. 

Motivation The two-channel Kondo impurity model 
consists of two identical species of non-interacting elec- 
trons antiferromagnetically coupled to a spin 1/2 impu- 
rity. Non-fermi liquid behavior results because of the 
inability to remove the impurity spin: it is energetically 
favorable for both conduction electron bands to couple 
to the impurity which gives a spin 1/2 complex on all 
length scales. As a result, the ground state is degenerate 
and the excitation spectrum non-Fermi liquid like. In 
contrast, the single channel Kondo model has a singlet 
ground state with the impurity spin screened out, and 
a Fermi liquid excitation spectrum corresponding to the 
removal of one conduction state from the system. On 
extension to the lattice and ignoring the renormalization 
of the environment around each spin, the array of single 
channel model singlets would simply renormalize the po- 
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tential scattering. In contrast, the array of many body 
doublets in the two-channel case would give rise to a dy- 
namical spin-disorder scattering in the absence of any 
cooperative transition that lent coherence to the spin ar- 
ray. We might thus anticipate a finite residual resistivity 
and other non-Fermi liquid properties in the paramag- 
netic phase 

Model The Hamiltonian for the two-channel Kondo 
lattice is 
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where c\ a (cj a a ) creates (destroys) an electron on site 
i in channel a = 1, 2 of spin a, Si is the Kondo spin on 
site z, and s.; jQ are the conduction electron spin operators 
for site i and channel a. The sites i form an infinite- 
dimensional hypercubic lattice. Hopping is limited to 
nearest neighbors with hopping integral t = t* j2\fd and 
the scaled hopping integral t* determines the energy unit 
and is set equal to one (t* = 1). Thus, on each site the 
Kondo spin mediates spin interaction between the two 
different channels. This problem is non trivial, and for 
the region of interest in which J > and T -C J, t* 
it is non perturbative. Clearly some simplifying method 
which allows for a solution of the lattice problem in a 
non-trivial limit is necessary. 

Formalism. Such a method was proposed by Metzner 
and Vollhardt Q who observed that the renormalizations 
due to local two-particle interactions become purely lo- 
cal as the coordination number of the lattice increases. 
More precisely, the irreducible single-particle self energy 
£g(z) and therefore the irreducible two-particle self en- 
ergy £ , -,(z, z') both become independent of momen- 
tum for large coordination number (2d — > oo) A 
further consequence of this is that the solution of most 
standard lattice models may be mapped onto the solu- 
tion of a local correlated system coupled to an effective 
bath that is self-consistently determined [TT|-[l^1. We re- 
fer the reader to the above references and recent reviews 
for further details on the method Gffl • 

Simulation. In order to solve the remaining impurity 
problem, we use the Kondo impurity algorithm of Fyc 
fyl , modified to simulate the two-channel problem jl9| . 
In the QMC the problem is cast into a discrete path 
formalism in imaginary time, where i~i — I At, At = 
(3/L, and L is the number of times slices. The values of 
L used ranged from 8 to 96, with the largest values of L 
reserved for the largest values of (3 since the time required 
by the algorithm increases like L? . A sign problem was 
also encountered in the QMC process and it also limited 
how low in temperature we could go. 

The Euclidean-time QMC results for the local greens 
function G(r) were then analytically continued to real 



frequencies using the "annealing" Maximum Entropy 
method [|lj. Here, we start from an exact result for 
the density of states (DOS) at high temperatures (i.e. 
a gaussian) and use this as a default model to analyti- 
cally continue the highest temperature data. The output 
is then used as a default model to analytically continue 
the data for the next lower temperature, and so on. As 
the temperature is lowered, only the low-energy features 
change, so very little additional information needs to be 
added at each temperature. To ensure that we have the 
best possible results from this procedure, we systemati- 
cally improve the precision of the QMC data (by running 
longer) until the continued spectra stabilize. The single 
particle self energy may then be obtained by inverting 
the relation G(oj) — —i^/nw (u> + fi — S(w)), where w(z) 
is the complex Fadeev function. 

Most of our simulations were limited to the model at 
half filling of the conduction band (N — 1.0 for J = 
0.75,0.625,0.5,0.4); however, we have also studied the 
system away from half filling (N — 0.75 and 0.50) for 
two values of J (J = 0.75 and 0.625). The qualitative 
features of the single-particle properties of the model do 
not depend strongly upon filling. 

Results. As the temperature is lowered, the local 
susceptibilities show evidence of Kondo screening: the 
screened local Kondo moment T\(T) falls like Tln(T), 
and the local Kondo spin-conduction spin correlation in- 
dicates an antiparallcl alignment (S^ • s^ Q ) < 0; whereas 
local spin-spin correlation function between the chan- 
nels indicates a weak parallel alignment (s^ • s^i) > 0. 
By comparing the susceptibility of the local moment to 
the impurity susceptibility, we are able to estimate the 
Kondo scale Tq from the defining relation Tqx(Tq) = 
0.3. We find that T Q ~ 0.85 Jexp(-1.01t*/ J). For 
comparison, in the impurity limit we would anticipate 
T™ p ~ 7r- 1 / 2 Jexp(-v^ : r/J). Hence, T is collectively 
enhanced in the lattice, as found in the one channel An- 
derson lattice problem near half filling P^-M ; we discuss 
our enhancement further below. 

When the temperature is lowered below T , we find 
non-fermi liquid behavior in the single-particle properties 
of the model. For example, in Fig. 1 the derivative of the 
particle distribution function n(tk) = T^2 n G(tk,ioJ n ) 
is plotted. This saturates to a finite width distribution, 
in contrast to a Fermi liquid where it would display a 
dominant delta function contribution as T — > 0, and a 
Luttinger liquid || or Marginal fermi liquid 20 1 where it 



would have a singular divergence. Clearly we see no such 
features in Fig. 1. 

In Fig. 2 we show one electron properties of the model. 
Fig. 2(a) displays the single particle DOS. This has a 
finite value as u>, T — > 0, with a peak away from the 
Fermi energy. Note that suppression of N(ui — 0) leads 
to an enhanced effective medium DOS which accounts 
qualtitatively for our collectively enhanced To values. 

Novel behavior is seen in the real part of the one 
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electron self energy which has positive slope at u — > 
(Fig. 2(c)). For a Fermi liquid, this slope would be 
negative. The physical content is important: Z = 
1/(1 — dReT,/du>) measures the overlap of the quasipar- 
ticle wave function with the original one-electron wave 
function having the same quantum labels. A positive 
slope leads to Z > 1 or Z < 0, indicating a breakdown 
of the quasiparticle concept. Concommitant with the fi- 
nite width of —dn/dek is a finite imaginary part to the 
low temperature self energy (Fig. 2(b)). This indicates 
that the one-electron excitations are ill defined on ap- 
proach to the Fermi surface, again ruling out a Fermi 
liquid description. Since the low temperature thermody- 
namic properties such as the specific heat, uniform mag- 
netic susceptibility, and charge susceptibility display no 
evidence for a gap, we believe the observed behavior in- 
dicate a new kind of metallic state. 

From our numerical results, one can show that the zero 
temperature self energy must be non-analytic. If the self 
energy is analytic everywhere on the real axis (and the 
upper half complex plane), then one may easily show that 
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where luq = ttT is the lowest fermionic Matsubara fre- 
quency along the imaginary axis. However, we find that 
when T < Tq, Im ^"° - < and is apparently divergent; 

whereas, from Fig. 3 it is apparent that dR ° s ^) > l. 

u— 

Thus, S(w) cannot be analytic at the origin of the com- 
plex plane. Given this, and the apparent shape of E(cj) 
in Fig. 2(b), we postulate the form ImS(w) ~ — c + \ui\. 

The non-fermi liquid behavior also strongly effects ex- 
perimentally relevant transport properties of the system. 
The electrical resistivity is shown in Fig. 3. We find that 
p(T)/p(0) curves for different J values collapse onto a 
universal scaling curve when plotted against T/Tq. As 
shown in the inset, p(T) ~ p(0)[l + B{T/T a )} for T -> 0, 
with B < |25[ . We interpret the finite value of p(0) to- 
gether with — im£(0, 0) as "spin disorder scattering" off 
of the degenerate screening clouds centered about each 
local moment spin. Given the finite density of one par- 
ticle excitations at the Fermi energy, this finite residual 
resistivity is indicative of an "incoherent metal phase" 
brought about by the disordered spin degrees of free- 
dom, in qualitative agreement with results obtained with 
a Lorentzian bare conduction DOS (which doesn't self- 
consistent ly renormalize) J8|. We conjecture that the an- 
tiferromagnetism at half filling will lift the disorder at 
each site and produce an insulating state due to the cell 
doubling. 

We have also computed the optical conductivity (not 
shown) and find that it displays vanishing Drude weight 
at low temperature together with a finite frequency peak. 
Both these features again support our interpretation in 
terms of a new kind of non-Fermi liquid metallic state. 



We find a commensurate antiferromagnetic transition 
at half filling for the simple nearest neighbor hopping 
model. We can suppress this phase by adding next neigh- 
bor hopping along each principle axis which will not alter 
our conclusions for the paramagnetic phase. The anti- 
ferromagnetic order is suppressed away from half filling. 
The complete phase diagram of the model will be ex- 
plored in a separate publication. 

Finally, we mention the possible applicability of our 
results to concentrated heavy fermion systems. Three 
systems display resistivity of the form p(T) w ,o(0)[l + 
B(T/T ) a ] for T < T , with p(0) of order the unitarity 
limit and a » 1. They are: ITCug^Pd* (with B < 0) [|, 
UBei 3 (B > 0) @, and CeCu 2 Si 2 (B > 0) ||. While 
UCus-^Pda; is a possible example of the distribution of 
Kondo scales scenario §, UBei 3 and CeCu 2 bi 2 are or- 
dered compounds which have been proposed as possible 
two-channel lattice systems (see Refs. @(b,c)); it is con- 
ceivable that 1 jd corrections can effect the sign change of 
B relative to our results. Another key difference with our 
model is that we assume global SU (2) \ sp in X SU (2) \ c hannei 
symmetry, while in the real materials these can only be 
local symmetries. 
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FIG. 1. Derivative of the particle distribution function 
n ( e *0 = TY, n G{€ k ,ioj n ) when J = 0.625 and N = 1.0. In a 
non-interacting system or one with a fermi-liquid low temper- 
ature state, —dn(0) / 'de^ oc Z 2 f3 and the distribution should 
have a width oc T which is the width of the fermi edge and is 
proportional to the scattering rate at the fermi surface. How- 
ever, in this case —dn(0)/de k saturates at low T and so does 
the width of —dn(ek)/dek , shown in the inset, indicating that 
a fermi liquid has not formed and that scattering rate remains 
finite when e fe = fi as T — > 0. 
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FIG. 2. Single particle properties of the two-channel Kondo 
lattice in infinite dimensions when J = 0.750 and N — 1.0. (a) 
Single-particle density of states (DOS). At high temperatures 
T ^> To (not shown), the DOS is a gaussian, crossing over to 
the peaked distribution with relative suppression at to — > 
for lower temperatures when T <g To. (b) Imaginary part of 
the self energy. As the temperature is lowered the self energy 
does not approach a fermi liquid form ImE(u) cx — T 2 — lo 2 , 
but rather appears to approach the non-analytic form (see 
text) ImS(w) cx — c + |w|. (c) The real part of the self en- 
ergy ReE(oj), is also anomalous since its initial slope is posi- 
tive indicating a quasiparticle renormalization factor which is 
greater than one. 
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FIG. 3. Resistivity of the two-channel Kondo lattice. 
p(T)/p(0) is plotted versus T/Tq for various values of J. In 
the inset, the lowest temperature data (for T/To < 1, shown 
as open circles) was fit to p(T)/p(0) = 1 + B (T/T ) a , with 
B = -0.4 and a = 1.03. 
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